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ABSTRACT 

A  series  of  flight  tests  were  conducted  on  the  PC  9/A  aircraft,  A23-045, 
at  the  Royal  Australian  Air  Force’s  Aircraft  Research  and  Development  Unit. 
System  identification  techniques  were  applied  to  the  data  obtained  from  these 
flight  tests  to  determine  the  stability  and  control  derivatives  of  the  aircraft. 
The  longitudinal  and  lateral  results  for  the  aircraft  in  approach  and  departure 
configurations  are  presented  in  this  report  and  comparisons  are  made  with 
wind  tunnel  estimates. 
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A  Correlation  between  Flight-determined  Longitudinal  and 
Lateral  Derivatives  and  Wind  Tunnel  Data  for  the  Pilatus 
PC  9/A  Training  Aircraft  in  Approach  and  Departure 

Configurations 


EXECUTIVE  SUMMARY 

Air  Operations  Division  (AOD)  has  developed,  or  acquired,  a  number  of  fixed-wing 
flight  dynamic  models  and  is  also  responsible  for  providing  advice  to  the  Australian  De¬ 
fence  Organisation  (ADO)  on  flight  simulator  flight  dynamic  model  requirements.  The 
models  generally  make  use  of  extensive  static  and  dynamic  stability  and  control  derivative 
databases,  in  addition  to  engine  and  flight  control  models.  The  static  model  data  may  be 
obtained  from  wind  tunnel  testing,  whilst  the  dynamic  data  is  traditionally  obtained  from 
flight  tests  using  system  identification  techniques. 

The  system  identification  techniques  used  to  estimate  aerodynamic  derivatives  of  con¬ 
ventional  aircraft  are  well  established.  The  major  requirement  of  these  techniques  is  high 
fidelity  measurements  of  manoeuvre  input  (i.e.  control  surface  deflections),  and  aircraft 
response  (i.e.  angular  rates  and  linear  accelerations),  as  well  as  air  data  including  airspeed, 
altitude,  angle-of-attack  and  angle-of-sideslip. 

Following  the  completion  of  the  AOD  PC  9/A  wind  tunnel  tests,  the  requirement 
existed  for  a  dedicated  flight  dynamic  modelling  flight  test  program  to  both  validate  the 
flight  dynamic  model  of  the  PC  9/ A  and  to  provide  dynamic  derivative  estimates  that 
were  not  obtainable  in  the  AMRL  wind  tunnel.  An  instrumentation  suite,  including 
an  air  data  probe  for  the  direct  and  accurate  measurement  of  angle-of-attack,  angle-of- 
sideslip,  temperature,  and  static  and  dynamic  pressure,  was  fitted  to  the  aircraft  and  a  test 
manoeuvre  matrix  was  designed  specifically  for  the  purpose  of  gathering  flight  dynamic 
data.  This  flight  test  program  was  conducted  at  the  Aircraft  Research  and  Development 
Unit  (ARDU)  during  1998/99. 

This  report  details  the  analysis  of  the  longitudinal  and  lateral  manoeuvres  carried  out 
with  the  aircraft  in  the  approach  and  departure  configurations.  The  static  and  dynamic 
derivatives  thus  obtained  are  compared  with  wind  tunnel  estimates.  A  discussion  of  some 
of  the  difficulties  encountered  during  the  estimation  process  is  also  included.  The  data 
obtained  from  these  tests  will  be  used  in  the  development  of  the  PC  9/ A  flight  dynamic 
model,  further  enhancing  AOD  support  for  the  PC  9/A  fleet,  including  possible  upgrades 
to  the  part-task  trainer. 


iii 


DSTO-TR-0479 


DSTO-TR-0479 


Authors 


Nick  van  Bronswijk 

Air  Operations  Division 

Nick  van  Bronswijk  graduated  from  the  University  of  Sydney 
with  first  class  honours  in  Aeronautical  Engineering  in  1995.  He 
is  currently  completing  his  PhD  in  the  area  of  propeller  power 
effects.  Since  joining  AMRL  in  1998  he  has  been  involved  with  a 
range  of  projects  including  flight  and  wind  tunnel  test  programs 
of  the  Pilatus  PC  9/ A. 


Andrew  Snowden 

Air  Operations  Division 

Andrew  Snowden  graduated  from  the  Royal  Melbourne  Insti¬ 
tute  of  Technology  in  1991,  having  obtained  a  Bachelor  of  Engi¬ 
neering  in  Aerospace  Engineering  with  first  class  honours.  Be¬ 
fore  commencing  work  at  AMRL  he  served  a  short  term  attach¬ 
ment  with  McDonnell  Douglas  in  St.  Louis,  Missouri,  as  the 
recipient  of  the  McDonnell  Aircraft  Company  Aerospace  Engi¬ 
neering  Prize  for  1991.  He  commenced  employment  at  AMRL 
in  1992  and  participated  in  several  experimental  investigations 
in  a  number  of  different  fields  during  his  rotation.  Now  a  full 
time  member  of  the  flight  dynamics  and  performance  area  of 
the  Air  Operations  Division,  his  most  recent  work  includes  the 
organisation  of  flight  testing  for  parameter  identification  of  the 
Pilatus  PC  9/ A  as  well  as  performance  estimation  of  modern 
jet  fighters. 


Hilary  Keating 

Air  Operations  Division 

Hilary  Keating  graduated  from  the  University  of  Sydney  in 
1998,  having  obtained  a  Bachelor  of  Engineering  in  Aeronauti¬ 
cal  Engineering  with  first  class  honours.  She  commenced  em¬ 
ployment  at  AMRL  in  1999,  and  has  been  involved  in  perfor¬ 
mance  modelling  and  parameter  identification  of  the  Pilatus 
PC  9/A. 


DSTO-TR-0479 


GeofF  Brian 

Air  Operations  Division 

Geoff  Brian  commenced  employment  at  the  Aeronautical  and 
Maritime  Research  Laboratory  in  1989,  after  graduating  from 
the  University  of  New  South  Wales  with  a  degree  in  Aeronauti¬ 
cal  Engineering,  Honours  Class  1.  During  his  employment  with 
AMRL,  he  has  been  involved  with  the  analysis  of  flight  data  for 
the  development  of  a  comprehensive  flight  dynamic  model  of  the 
F-lllC  aircraft;  the  development  of  aircraft  performance  esti¬ 
mation  applications;  and  involvement  with  a  performance  flight 
trial  of  a  F-lllC  aircraft  fitted  with  Pratt  &  Whitney  TF30-P- 
109  engines.  Mr.  Brian  was  attached  to  the  Royal  Australian 
Air  Force  -  Aircraft  Research  and  Development  Unit  in  1992, 
where  he  was  involved  with  a  performance  flight  trial  of  an 
air-to-air  refuelling  capable  Boeing  707  aircraft.  During  1998 
and  1999  Mr.  Brian  attended  Loughborough  University  in  the 
United  Kingdom  being  awarded  with  a  MSc  in  Industrial  Math¬ 
ematical  Modelling  with  Distinction.  Mr  Brian  is  currently  the 
manager  for  aircraft  flight  dynamic  activities  at  AMRL. 


vi 


DSTO-TR-0479 


Contents 

Glossary  xi 

Notation  xi 

1  Introduction  1 

2  PC  9/ A  Test  Aircraft  1 

2.1  Aircraft  Description .  1 

2.2  Flight  Control  System .  2 

2.3  Weight,  Centre-of-Gravity  and  Mass  Moments-of-Inertia .  2 

3  Instrumentation  2 

4  Methods  of  Analysis  3 

4.1  Stepwise  Regression .  3 

4.1.1  Error  Band  .  4 

4.2  Maximum  Likelihood  .  4 

4.2.1  Cramer- Rao  Bounds .  6 

5  Results  and  Discussion  6 

5.1  Angle-of- Attack  Derivatives .  6 

5.2  Pitch  Rate  Derivatives .  7 

5.3  Longitudinal  Control  Derivatives .  7 

5.4  Angle-of-Sideslip  Derivatives  .  8 

5.5  Yaw  Rate  Derivatives .  8 

5.6  Roll  Rate  Derivatives  .  9 

5.7  Lateral  Control  Derivatives .  9 

5.8  Variation  of  Derivatives  with  Airbrake  Deployment .  9 

5.9  Configuration  and  Altitude  Dependency .  11 

6  Conclusions  13 

15 


vii 


References 


Appendices 


DSTO-TR-0479 


A  Mass,  Cent re-of- Gravity  and  Mass  Moments-of-Inertia  37 

Figures 

1  Air  data  boom  installation .  17 

2  PC  9/A  angle-of-attack  and  pitch  rate  derivatives,  gear  down .  18 

3  PC  9/A  angle-of-attack  and  pitch  rate  derivatives,  take-off  flap .  19 

4  PC  9/ A  angle-of-attack  and  pitch  rate  derivatives,  landing  flap . 20 

5  PC  9/A  longitudinal  control  derivatives,  gear  down . .  •  21 

6  PC  9/A  longitudinal  control  derivatives,  take-off  flap . 22 

7  PC  9/ A  longitudinal  control  derivatives,  landing  flap . 23 

8  PC  9/ A  angle-of-sideslip  derivatives,  gear  down .  24 

9  PC  9/ A  angle-of-sideslip  derivatives,  take-off  flap . 25 

10  PC  9/A  angle-of-sideslip  derivatives,  landing  flap . 26 

11  PC  9/ A  yaw  rate  derivatives,  gear  down . 27 

12  PC  9/ A  yaw  rate  derivatives,  take-off  flap . 28 

13  PC  9/ A  yaw  rate  derivatives,  landing  flap . 29 

14  PC  9/ A  roll  rate  derivatives,  gear  down .  30 

15  PC  9/ A  roll  rate  derivatives,  take-off  flap.  .  . .  31 

16  PC  9/A  roll  rate  derivatives,  landing  flap . 32 

17  PC  9/ A  lateral  control  derivatives,  gear  down .  33 

18  PC  9/ A  lateral  control  derivatives,  take-off  flap . 34 

19  PC  9/A  lateral  control  derivatives,  landing  flap . 35 

Tables 

1  Longitudinal  derivatives  estimated .  5 

2  Lateral  derivatives  estimated .  5 

3  Comparison  of  averaged  longitudinal  flight  test  derivatives  with  airbrake  de¬ 

ployed,  using  maximum  likelihood  and  stepwise  regression  techniques,  to  av¬ 
eraged  clean  aircraft  results  [1] .  10 

4  Comparison  of  averaged  lateral  flight  test  derivatives  with  airbrake  deployed, 

using  maximum  likelihood  and  stepwise  regression  techniques,  to  averaged 
clean  aircraft  results  [2] .  10 

5  Signal  variance  ratios  for  gear  down,  take-off  flap,  and  landing  flap  configu¬ 
rations .  11 

viii 


DSTO-TR-0479 


/ 


6  Summary  of  longitudinal  derivatives  estimated  using  maximum  likelihood.  .  .  12 

7  Summary  of  longitudinal  derivatives  estimated  by  stepwise  regression . 12 

8  Summary  of  lateral  derivatives  estimated  using  maximum  likelihood . 13 

9  Summary  of  lateral  derivatives  estimated  by  stepwise  regression . 13 

A1  Flight  test  aircraft  mass  distribution  [3,4] . 38 

A2  Retracted  landing  gear  mass  and  moment  arm . 38 

A3  Deployed  landing  gear  mass  and  moment  arm . 38 

A4  Mass  moment-of- inertia  increments  due  to  landing  gear  deployment . 39 


IX 


DSTO-TR-0479 


X 


Glossary 


AHRS  Artificial  Horizon  Reference  System 

AMRL  Aeronautical  and  Maritime  Research  Laboratory 

AOD  Air  Operations  Division 

ARDU  Aircraft  Research  and  Development  Unit 

DSTO  Defence  Science  and  Technology  Organisation 

GDAS  General  Data  Acquisition  System 

IAS  Indicated  Airspeed 

ITT  Inlet  Turbine  Temperature 

KIAS  Knots  Indicated  Airspeed 

MAC  Mean  Aerodynamic  Chord 

NASA  National  Aeronautics  and  Space  Administration 

NG  Gas  generator  speed 

OAT  Outside  Air  Temperature 

PAP  Primary  Analysis  Processor 

RAAF  Royal  Australian  Air  Force 

SHP  Shaft  Horse  Power 


DSTO-TR-0479 


xii 


Notation 

ax,0’Yi<^z  Body  axes  linear  accelerations  (g) 

b  Reference  span  (10.124  m) 

Cl  Rolling  moment  coefficient 

Ci-W 

Ci^  Rolling  moment  coefficient  due  to  angle-of-sideslip 

degree) 

Qp  Rolling  moment  coefficient  due  to  roll  rate 

radian) 

Ci^  Rolling  moment  coefficient  due  to  yaw  rate 

(per  radian) 

Ci^  Rolling  moment  coefficient  due  to  aileron  deflection 

=  S  (P®'' 

Cig  Rolling  moment  coefficient  due  to  rudder  deflection 

^isr  =  (P®^  ‘^®g^®®)  _ 

Cm  Pitching  moment  coefficient 

—  qSc 

Cmo,  Pitching  moment  coefficient  due  to  angle-of-attack 

Cma  =  ^  (per  degree) 

Cma  Pitching  moment  coefficient  due  to  rate  of  change  of  angle-of-attack 

^rna  =  (per  radian) 

Cms  Pitching  moment  coefficient  due  to  elevator  deflection 

^  (P®’^  degree) 

Cmg  Pitching  moment  coefficient  due  to  pitch  rate 

Cmg  =  (per  radian) 

Cn  Normal  force  coefficient 

qS 

Cxa  Normal  force  coefficient  due  to  angle-of-attack 

^  (per  degree) 

Cn^  Normal  force  coefficient  due  to  rate  of  change  of  angle-of-attack 

C'Afd  =  (per  radian) 

Cjv^  Normal  force  coefficient  due  to  elevator  deflection 

w  (per  degree) 

Cjv  Normal  force  coefficient  due  to  pitch  rate 

C'n,  =  (per  radian) 

Cn  Yawing  moment  coefficient 

_  n 

'-'n  -  qSb 

Cnp  Yawing  moment  coefficient  due  to  angle-of-sideslip 

^  (per  degree) 

Cnp  Yawing  moment  coefficient  due  to  roll  rate 

Cnp  =  (per  radian) 

Cur  Yawing  moment  coefficient  due  to  yaw  rate 

^nr  =  (per  radian) 


DSTO-TR-0479 


Cn^  Yawing  moment  coefficient  due  to  aileron  deflection 

Cns  Yawing  moment  coefficient  due  to  rudder  deflection 

^  W  (P^^  degree) 

Cx  Longitudinal  force  coefficient 

Cx  =  § 

Cy  Side  force  coefficient 

Cy^  Side  force  coefficient  due  to  angle-of-sideslip 

^  (P®^  degree) 

Cyp  Side  force  coefficient  due  to  roll  rate 

==  (P®^  radian) 

Cy^  Side  force  coefficient  due  to  yaw  rate 

^Yr  =  (per  radian) 

Cy^  Side  force  coefficient  due  to  aileron  deflection 

^Y,,  =  ^  (per  degree) 

Cy^^  Side  force  coefficient  due  to  rudder  deflection 
C'n,  =  ®  (per  degree) 

Cz  Vertical  force  coefficient 

Cz  =  § 

c  Reference  chord  (1.65  m) 

F  Ratio  of  the  regression  mean  square  to  the  residual  mean  square 
g  Gravitational  acceleration  (9.81  m/s^) 

j^xx  Roll  moment  of  inertia  (kg.m^) 

Jyy  Pitch  moment  of  inertia  (kg.m^) 
j^zz  Yaw  moment  of  inertia  (kg.m^) 

IxY  Cross  product  of  moment  of  inertia  (kg.m^) 

Ixz  Cross  product  of  moment  of  inertia  (kg.m^) 

lyz  Cross  product  of  moment  of  inertia  (kg.m^) 

I  Rolling  moment  (N.m) 

M  Mass  of  aircraft  (kg) 

m  Pitching  moment  (N.m) 

N  Normal  force  (N) 

n  Yawing  moment  (N.m) 

p  Roll  rate  (rad/s) 

q  Pitch  rate  (rad/s) 

q  Dynamic  pressure  (N/m^) 

P?  Squared  multiple  correlation  coefficient 
r  Yaw  rate  (rad/s) 

S  Reference  area  (16.29  m^) 

V  True  velocity  (m/s) 

X  Longitudinal  force  (N) 

Xr  a  Longitudinal  c.q.  position  (m) 

Y  Side  force  (N) 

Uc.g.  Lateral  c.g.  position  (m) 

Z  Vertical  force  (N) 

Zc.g.  Vertical  c.g.  position  (m) 

xiii 


DSTO-TR-0479 


a  Angle-of-attack  (°) 

a  Time  derivative  of  angle-of-attack  (rad/ s) 

I3  Angle-of-sideslip  (°) 

5a  Aileron  deflection  =  5a^  —  {°) 

5a  Elevator  deflection  (°) 

5r  Rudder  deflection  (°) 

0  Pitch  angle  (°) 

(p  Roll  angle  (°) 

Ip  Yaw  angle  (°) 

Statistical  variance 

Subscripts 

b  Body  axes 
c.g.  Centre-of-gravity 
dep  Deployed 
F  Fuel 
GD  Gear  down 
L  Left  (port) 
mg  Main  gear 
ng  Nose  gear 
ret  Retracted 
R  Right  (starboard) 
tot  Total 


XIV 


DSTO-TR-0479 


Aircraft  Principle  Dimensions. 

(Numeric  data  sourced  from  Pilahis  Stnictural  Configuration  Drawing  506.00.09.220F) 


DSTO-TR-0479 


1  Introduction 

The  Pilatus  PC  9/A  is  one  of  a  number  of  high  performance  turbo-prop  aircraft  currently 
operated  by  the  Royal  Australian  Air  Force  (RAAF).  Work  currently  underway  in  the  Air 
Operations  Division  (AOD)  of  the  Defence  Science  and  Technology  Organisation  (DSTO) 
includes  investigations  into  the  propeller  power  effects  of  such  aircraft,  and  the  PC  9/A 
was  considered  a  suitable  platform  for  study.  A  six  degree-of- freedom  flight  dynamic  model 
of  the  PC  9/A  has  been  developed  by  AOD  for  use  in  these  investigations,  as  well  as  for 
pilot-in-the-loop  simulations  in  the  Air  Operations  Simulation  Centre  and  incident  and 
accident  investigations. 

The  development  of  a  flight  dynamic  model  requires  information  on  the  static  and  dynamic 
stability  and  control  derivatives  of  the  aircraft,  as  well  as  flight  control  laws  and  physical 
properties.  The  static  data  for  the  PC  9/A  were  collected  during  both  power-off  and  power- 
on  testing  of  a  scaled  aircraft  model  in  the  Aeronautical  and  Maritime  Research  Laboratory 
(AMRL)  low-speed  wind  tunnel.  Additional  power-off  static  data  were  obtained  from  a 
computational  fluid  dynamic  model  for  a  limited  number  of  aircraft  configurations. 

Flight  test  data  for  the  PC  9/ A  was  required  to  both  validate  the  fiight  dynamic  model 
and  to  provide  dynamic  derivative  estimates  of  the  aircraft.  A  fiight  test  program  was  con¬ 
ducted  on  a  PC  9/ A  aircraft,  serial  number  A23-045,  by  the  RAAF  Aircraft  Research  and 
Development  Unit  (ARDU)  between  November  1998  and  February  1999.  The  details  of 
this  fiight  test  program  and  the  aircraft  instrumentation  are  reported  in  reference  [5].  The 
estimation  of  the  stability  and  control  derivatives  from  the  flight  test  data  involved  the  use 
of  system  identification  techniques,  and  required  control  input,  aircraft  response  and  flight 
condition  data,  measured  using  a  high  fidelity  instrumentation  system.  Maximum  likeli¬ 
hood  and  stepwise  regression  techniques  were  employed  for  the  system  identification.  The 
longitudinal  derivatives  for  the  PC  9/A  aircraft  in  cruise  configuration  determined  using 
system  identification  techniques  are  reported  in  reference  [1],  while  the  lateral  derivatives 
for  the  aircraft  in  cruise  configuration  are  reported  in  reference  [2]. 

This  report  presents  the  longitudinal  and  lateral  stability  and  control  derivatives  of  the 
PC  9/A  extracted  from  flight  test  data  for  the  approach  and  departure  configurations, 
including  a  comparison  with  AMRL  wind  tunnel  test  data.  Sections  2  and  3  present 
details  of  the  test  aircraft  and  instrumentation  system.  Section  4  presents  the  system 
identification  techniques,  while  section  5  discusses  the  results. 

2  PC  9/ A  Test  Aircraft 

2.1  Aircraft  Description 

The  flight  test  aircraft,  A23-045,  is  a  Pilatus  PC  9/A  operated  by  ARDU.  The  PC  9/A  is 
a  single-engine,  metal-skinned,  low-wing,  tandem  two-seat  training  aircraft.  The  aircraft 
is  powered  by  a  Pratt  and  Whitney  PT6A-62  turbo-prop  engine  flat  rated  to  950  SHP  [3], 
which  drives  a  Hartzell  HC-D4N-2A  four-blade  variable  pitch  propeller.  The  aircraft  was 
instrumented  as  outlined  in  section  3. 
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2.2  Flight  Control  System 

The  aircraft  primary  flight  controls  consist  of  the  ailerons,  rudder  and  elevator.  The  control 
surfaces  are  manually  operated  from  a  conventional  dual  control  column  and  rudder  pedal 
arrangement.  The  stick  and  rudder  pedals  are  connected  to  the  control  surfaces  through 
a  system  of  control  rods,  bellcranks,  cables  and  levers.  Trimming  control  is  provided  on 
all  three  axes.  The  aircraft  is  fitted  with  trailing  edge  split  flaps  that  cover  the  final  18% 
of  the  wing  chord  from  the  wing  root  to  the  inboard  edge  of  the  aileron.  These  flaps  have 
three  discrete  positions,  retracted  {5f  —  0°),  take-off  {6f  =  23°),  and  landing  {6f  =  50  ). 
Additionally,  the  aircraft  is  fitted  with  a  perforated  flat  plate  airbrake  that  is  located  under 
the  centre  fuselage  between  the  wing  flaps.  This  airbrake  is  either  retracted  or  deployed 
at  an  angle  of  50° . 


2.3  Weight,  Centre-of-Gravity  and  Mass  Moments-of-Inertia 

The  test  aircraft,  A23-045,  was  weighed,  in  clean  configuration,  by  ARDU  prior  to  the 
flight  tests  and  had  a  basic  mass  of  1784.5  kg  and  a  longitudinal  centre-of-gravity  po¬ 
sition  of  26.25%  Mean  Aerodynamic  Chord  (MAC)  [4].  During  the  flight  test  program, 
the  aircraft  weight,  centre-of-gravity  and  mass  moments-of-inertia  varied  with  fuel  us¬ 
age.  Equations  for  the  mass,  centre-of-gravity,  and  mass  moments-of-inertia  of  the  clean 
aircraft  are  documented  in  [1,  Appendix  A].  These  equations  were  used  for  analysis  of 
flight  test  cases  with  airbrake  deployed.  Appendix  A  develops  equations  for  the  mass, 
centre-of-gravity,  and  mass  moments-of-inertia  of  the  aircraft  in  approach  and  departure 
configurations. 


3  Instrumentation 

Aircraft  A23-045  was  fitted  with  an  instrumentation  system  designed  specifically  for  the 
gathering  of  flight  dynamic  data.  A  summary  of  the  instrumentation  used  in  the  flight 
test  program  is  included  below  and  a  more  comprehensive  description  of  the  design  re¬ 
quirements  and  calibration  is  included  in  references  [5]  and  [6].  The  ARDU  General  Data 
Acquisition  System  (GDAS)  was  fitted  in  place  of  the  rear  ejection  seat.  The  data  were 
encoded  by  a  16-bit  pulse  code  modulation  system  and  were  recorded  onboard  the  air¬ 
craft  using  a  MARS-2000  14-track  tape.  Real-time  flight  test  monitoring  was  provided  by 
telemetry  data  transmitted  to  the  ARDU  Primary  Analysis  Processor  (PAP)  hut. 

The  angular  rates  (p,  g,r)  and  linear  accelerations  {ax^ay^az)  were  measured  using  the 
ARDU  KAISG1134-1  Motion  Platform.  This  consisted  of  three  Smith  Industries  950  RGS 
angular  rate  gyros  and  three  SunStrand  QA1400  servo  accelerometers.  The  aircraft  roll, 
pitch  and  yaw  attitude  angles  {(f>,  9,  ip)  were  obtained  by  tapping  output  from  the  existing 
LISA  2000A  Artificial  Horizon  Reference  System  (AHRS). 

Static  outside  air  temperature  (OAT),  indicated  airspeed  (IAS),  angle-of-attack  (a)  and 
angle-of-sideslip  (^0)  were  obtained  from  the  Rosemount  Model  92AN  flight  test  air  data 
boom,  mounted  on  the  outboard  hardpoint  on  the  starboard  wing  (see  figure  1).  Aileron 
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and  elevator  deflections  were  measured  using  Space  Age  Control  Inc.  series  160  cable 
position  transducers.  Rudder  deflection  was  measured  using  a  type  26V-11CX4C  position 
transducer.  All  sensors  were  calibrated  prior  to  commencement  of  the  flight  test  program. 

The  aircraft  was  also  instrumented  to  measure  engine  torque,  propeller  speed,  inlet  turbine 
temperature  (ITT),  gas  generator  speed  (NG),  fuel  flow  and  fuel  quantity.  The  torque, 
propeller  speed,  altitude  and  true  airspeed,  when  used  in  conjunction  with  a  performance 
map  for  the  Hartzell  HC-D4N-2A  propeller  [7],  allowed  the  calculation  of  engine  thrust. 


4  Methods  of  Analysis 


During  the  approach  and  departure  phase  of  flight  test  program,  doublet  manoeuvres  were 
performed  about  steady  flight  conditions  at  airspeeds  between  120  and  200  KIAS.  Refer¬ 
ence  [5]  describes  these  manoeuvres  in  detail.  The  25  longitudinal  and  27  lateral  doublet 
manoeuvres,  for  the  aircraft  in  approach  or  departure  configurations,  were  analysed  using 
stepwise  regression  [8]  and  maximum  likelihood  estimation  techniques  [9]  to  determine  the 
longitudinal  and  lateral  stability  and  control  derivatives  of  the  aircraft. 

Aircraft  derivatives  estimated  from  AMRL  power-off  wind  tunnel  tests  [10]  were  used  to 
provide  comparisons  for  the  static  aircraft  derivatives  and,  in  the  case  of  the  maximum 
likelihood  analysis,  to  provide  a  priori  estimates  to  increase  the  rate  of  convergence  of  the 
system  identification.  Due  to  the  lack  of  dynamic  derivative  estimates  for  the  aircraft  in 
approach  or  departure  configurations,  empirical  estimates  of  the  clean  aircraft  dynamic 
derivatives  were  used  as  a  priori  values  to  improve  convergence. 

A  right  handed  orthogonal  axes  system  was  adopted  for  the  analysis  of  the  flight  test  data. 
Positive  control  surface  deflections  were  defined  as  elevator  trailing  edge  down,  rudder 
trailing  edge  to  port,  and  starboard  aileron  trailing  edge  down,  port  aileron  trailing  edge 
up. 

4.1  Stepwise  Regression 

Stepwise  regression  is  an  unbiased  least  squares  estimator  in  which  new  independent  vari¬ 
ables  are  inserted  into  a  model,  one  at  a  time,  until  the  regression  equation  is  deemed 
acceptable.  The  appropriateness  of  the  model  can  be  determined  by  examining  a  num¬ 
ber  of  quantities  including  the  squared  multiple  correlation  coefficient  and  the  F  statistic. 
The  squared  multiple  correlation  coefficient,  B?,  gives  a  measure  of  the  importance  of 
each  variable  as  it  is  inserted  into  the  equation  [8];  however,  the  improvement  in  due 
to  the  addition  of  new  terms  must  have  some  real  significance  besides  simply  reflecting  the 
inclusion  of  more  terms.  This  can  be  determined  by  monitoring  the  F  statistic,  the  ratio  of 
the  regression  mean  square  to  the  residual  mean  square.  The  inclusion  of  any  significant 
terms  is  generally  accompanied  by  an  increase  in  the  F  statistic  and  the  best  fit  with  the 
least  number  of  parameters  may  be  obtained  by  maximising  F.  Any  variable  which  does 
not  make  a  significant  contribution  is  removed  from  the  model,  with  the  selection  process 
continuing  until  no  new  variables  remain  to  be  inserted  into  the  equation.  Whilst  stepwise 
regression  gives  estimates  of  the  derivatives  included  in  the  regression  equation,  it  also 
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permits  a  suitable  structure  for  the  model  to  be  determined.  In  addition,  it  provides  an 
independent  check  on  the  data  estimated  using  the  maximum  likelihood  technique. 

The  stepwise  regression  technique  was  applied  using  code  available  in  the  MATLAB© 
Statistics  Toolbox  [11].  The  following  model  equations  were  identified  during  the  analysis. 


Cn  —  Cno  +  (1) 

Cm  =  Cmo  +  Cma<^  +  ^  +  Cms^  (2) 

Cy  =  Cyq  +  Cy^^  +  Cy^^  5r  (3) 

C'n  =  Cno  +  Cnp^  +  C'np  ^  ^  +  Cng^  <^r  (4) 

Cl  =  Cl,  +  Q,/3  +  Ci^^  +  Ci^—  +  Ci,Ja  (5) 


4.1.1  Error  Band 

Included  in  the  stepwise  regression  analysis  is  the  calculation  of  an  error  band  on  the 
estimated  derivatives.  For  a  confidence  interval  of  95%,  this  error  band  is  approximately 
equal  to  two  standard  deviations.  Figures  2  to  19  show  the  stepwise  regression  derivative 
estimates,  including  the  calculated  error  band. 


4.2  Maximum  Likelihood 

The  maximum  likelihood  estimation  technique  was  applied  using  the  computer  program, 
pEst,  developed  at  NASA  Dryden  Flight  Research  Center  [12].  pEst  is  an  interactive  pa¬ 
rameter  estimation  program  which  solves  a  vector  set  of  time- varying,  ordinary  differential 
equations  of  motion. 

The  longitudinal  derivatives  identified  using  the  maximum  likelihood  technique  are  given 
in  table  1.  Three  aerodynamic  derivatives  used  in  the  longitudinal  flight  dynamic  model 
were  not  identified  due  to  difficulties  in  estimation  and  their  small  contribution. 

The  derivative  Cma  difficult  to  estimate  from  standard  flight  test  manoeuvres  and  is 
usually  considered  as  an  additional  component  of  the  pitch  damping  derivative  Cmg  during 
parameter  identification  [13].  Experience  gained  during  the  analysis  of  F-lllC  flight  test 
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Normal  Force 

Pitching  Moment 

Aerodynamic 

Cno. 

CrUa 

Cn.. 

Control 

Table  1:  Longitudinal  derivatives  estimated. 

data  [14]  has  shown  that  it  is  preferable  to  fix  Cm^  at  it’s  a  priori  value,  chosen  here  as 
-7.823  (per  radian),  and  allow  Cmg  to  be  estimated. 

Difficulties  also  arise  when  trying  to  estimate  both  the  normal  force  due  to  pitch  rate 
derivative,  C^v^,  and  the  normal  force  due  to  the  time  rate-of-change  of  angle-of-attack 
derivative,  Cn^.  Previous  analysis  of  the  PC  9/A  clean  aircraft  longitudinal  data  [1, 
Appendix  B]  has  shown  that  the  relative  contributions  of  and  to  the  total  aircraft 
normal  force  are  negligible.  Therefore,  these  derivatives  have  been  fixed  at  their  a  priori 
values,  chosen  as  7.96  (per  radian)  and  3.061  (per  radian),  respectively,  without  any  loss 
of  accuracy  of  the  estimation  of  the  other  normal  force  derivatives,  Cnc  and  . 

Due  to  a  lack  of  empirical  data  for  the  PC  9/ A  in  approach  and  departure  configurations, 
the  a  priori  values  of  ^Ng  and  Cat^,  used  during  the  data  analysis,  are  the  clean 

aircraft  values  obtained  from  reference  [15].  The  resulting  total  force  and  moment  coef¬ 
ficient  equations  used  in  the  longitudinal  state  and  response  equations  for  the  maximum 
likelihood  analysis  of  the  PC  9/ A  are  given  in  equations  6  and  7. 


no  cxc 

Cm  =  Cmo  +  Cmu^  +  +  3.061—  (6) 


no  (xc 

Cm  =  Cmo  +  Cma  “  +  C'm,  ^  +  Cmg^  Sg  -  7.823—  (7) 

The  lateral  derivatives  identified  using  the  maximum  likelihood  technique  are  given  in 
table  2. 


Side  Force 

Yawing  Moment 

Rolling  Moment 

Aerodynamic 

Cyu 

Cn, 

Cl, 

Cn, 

Cl, 

Cnr 

Cl. 

Control 

Cn,, 

Ci,„ 

Table  2:  Lateral  derivatives  estimated. 

Five  aerodynamic  derivatives  used  in  the  lateral  flight  dynamic  model  were  not  identified 
due  to  both  their  small  contribution  and  difficulty  of  estimation.  Cyj,  was  fixed  at  a  value 
of -0.16  (per  radian),  Cvr  at  0.34  (per  radian),  Cys^  at  0.0  (per  degree),  C^^  at  -0.000035 
(per  degree)  and  Ci^^  remained  fixed  at  0.0035  (per  degree).  As  for  the  longitudinal  case, 


5 


DSTO-TR-0479 


a  lack  of  approach  and  departure  configuration  empirical  data  required  the  use  of  clean 
aircraft  a  priori  values  for  Cyp,  Cy^,  obtained  from  references  [16] 

and  [15].  The  resulting  total  force  and  moment  coefficient  equations  used  in  the  lateral 
state  and  response  equations  for  the  maximum  likelihood  analysis  of  the  PC  9/A  are  given 
in  equations  8,  9  and  10.  The  diflFerences  between  equations  1,  2,  3,  4  and  5,  and  6,  7,  8,  9 
and  10  are  important,  and  arise  from  the  inclusion  of  the  a  priori  values  in  the  maximum 
likelihood  technique. 


Tyb  T'b 

Cy  =  Cy.  +  Cy„/3  -  0.1608^  +  0.3417^  +  Cy,/.  (8) 

c„  =  C„„  +  +  C„, ^  +  Cn,^-  0.0000353.  +  3.  (9) 

C,  =  C,.  +  a,/3  +  Q,  ^  +  C,.  ^  +  a.  3.  +  0.0034913,  (10) 


4.2.1  Cramer-Rao  Bounds 

For  the  estimated  parameters,  pEst  calculates  a  measure  of  the  estimation  certainty  known 
as  the  Cramer-Rao  bound.  A  detailed  interpretation  of  this  quantity  is  given  in  [17].  The 
Cramer-Rao  bounds  are  shown  for  each  derivative  estimated  by  pEst  in  figures  2  to  19. 
The  Cramer-Rao  bounds  have  been  factored  in  accordance  with  the  procedures  described 
in  [17]  to  account  for  the  presence  of  band-limited  noise. 


5  Results  and  Discussion 

Longitudinal  and  lateral  derivatives  estimated  via  maximum  likelihood  parameter  esti¬ 
mation,  stepwise  regression,  and  power-off  wind  tunnel  experiments  are  summarised  in 
figures  2  to  19.  Data  are  presented  for  the  aircraft  with  gear  down,  take-off  flap  with  gear 
down,  landing  flap  with  gear  down,  and  clean  with  airbrake  extended. 


5.1  Angle-of- Attack  Derivatives 

Figures  2  to  4  show  the  angle-of-attack  derivatives  and  Cma  plotted  against  a.  The 
flight  test  estimates  of  are  generally  in  good  agreement  with  wind  tunnel  estimates. 
The  maximum  likelihood  estimates  are  higher  than  the  stepwise  regression  estimates  and 
there  is  an  increase  in  the  Cramer-Rao  bounds  of  the  estimates  as  flap  deflection  is  in¬ 
creased,  indicating  increased  uncertainty.  This  is  discussed  in  more  detail  in  section  5.9. 

The  flight  test  estimates  of  pitch  stiffness.  Cm,,,  are  consistently  smaller  in  magnitude  than 
the  wind  tunnel  estimates,  which  are  power-off.  This  difference  is  due  to  the  imbalance 
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between  two  competing  effects  acting  on  the  portion  of  the  horizontal  tail  within  the 
propeller  slipstream,  these  being  the  increased  dynamic  pressure  and  the  reduced  angle  of 
incidence.  The  reduction  in  angle  of  incidence  dominates,  resulting  in  an  overall  reduction 
in  the  pitch  stiffness  of  the  aircraft.  This  phenomenon  is  discussed  in  more  detail  in 
reference  [1].  The  stepwise  regression  estimates  are  consistently  smaller  in  magnitude 
than  the  maximum  likelihood  estimates.  The  small  Cramer-Rao  bounds  indicate  a  high 
level  of  confidence  in  the  maximum  likelihood  estimates. 


5.2  Pitch  Rate  Derivatives 

Figures  2  to  4  show  the  pitch  rate  derivative,  Cmg,  plotted  against  a.  As  discussed  earlier, 
Cm^  was  fixed  at  it’s  a  priori  value  to  aid  the  estimation  of  the  pitch  damping  derivative 
Crag  ■  This  method  has  been  used  successfully  during  other  parameter  estimation  exercis¬ 
es  [14],  and  was  verified  as  suitable  for  use  with  the  PC  9/A  flight  testing  in  reference  [1]. 

The  AMRL  low-speed  wind  tunnel  has  no  facility  to  measure  dynamic  derivatives,  and 
empirical  data  was  not  available  for  the  Pilatus  PC  9/ A  in  approach  or  departure  con¬ 
figurations.  Figures  2  to  4  show  that  the  maximum  likelihood  estimates  are  smaller  in 
magnitude  than  the  stepwise  regression  estimates,  and  that  the  Cramer-Rao  bounds  and 
difference  between  maximum  likelihood  and  stepwise  regression  estimates  increases  as  flap 
deflection  increases. 


5.3  Longitudinal  Control  Derivatives 

Figures  5  to  7  show  the  longitudinal  control  derivatives,  Cms^  and  ,  plotted  against  a. 
The  maximum  likelihood  estimates  of  Cms^  agree  very  closely  with  wind  tunnel  estimates 
for  all  configurations.  The  stepwise  regression  estimates  are  smaller  in  magnitude.  Small 
Cramer-Rao  bounds  indicate  high  confidence  in  the  estimates,  although  they  increase  with 
increasing  flap  deflection. 

Estimates  of  C^s  show  only  a  small  amount  of  scatter  for  all  configurations.  Maximum 
likelihood  estimates  are  positive,  and  just  below  the  wind  tunnel  estimates,  while  stepwise 
regression  estimates  are  negative.  As  discussed  previously,  a  priori  values  are  included 
in  the  maximum  likelihood  technique,  whereas  the  stepwise  regression  technique  does 
not  consider  these  values.  Analysis  of  the  effect  of  inclusion  or  exclusion  of  the  a  priori 
values  [1]  demonstrated  a  change  in  sign  of  the  maximum  likelihood  estimates  of  .  It 
was  therefore  theorised  that  the  effect  of  setting  both  and  Cyv,  to  zero  is  to  reduce  the 
rate  of  onset  of  normal  force  during  a  manoeuvre  and  therefore,  the  modelling  techniques 
decrease  the  value  of  ,  to  below  zero,  to  compensate.  This  ensures  the  overall  match  of 
the  data  is  good  but  results  in  some  loss  of  model  resolution  during  the  initial  application 
of  elevator  [1]. 
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5.4  Angle-of-Sideslip  Derivatives 

Figures  8  to  10  show  the  angle-of-sideslip  derivatives  Cn^  and  Ci^  plotted  against  a. 
The  flight  test  estimates  of  Cy^  are  greater  in  magnitude  than  the  wind  tunnel  estimates, 
and  the  difference  between  flight  test  and  wind  tunnel  increases  with  flap  deflection.  The 
Cramer-Rao  bounds  also  increase  with  flap  deflection.  In  general,  the  stepwise  regression 
estimates  are  smaller  in  magnitude  than  the  maximum  likelihood  estimates. 

Flight  test  estimates  of  the  directional  stability  derivative,  Cnp ,  show  good  agreement  be¬ 
tween  the  stepwise  regression  and  maximum  likelihood  techniques  for  the  gear  down  cases, 
but  a  slight  divergence  between  the  two  techniques  as  flap  deflection  is  increased.  Both 
flight  test  estimation  techniques  yield  values  significantly  below  that  measured  in  the  wind 
tunnel,  which  is  power-off.  This  difference  is  due  to  the  imbalance  between  two  competing 
effects  acting  on  the  portion  of  the  vertical  tail  within  the  propeller  slipstream,  these  being 
the  increased  dynamic  pressure  and  the  decreased  angle  of  incidence.  The  reduction  in 
angle  of  incidence  dominates,  resulting  in  an  overall  reduction  in  the  directional  stability 
of  the  aircraft. 

The  flight  test  estimates  of  the  dutch  roll  derivative,  show  good  agreement  with 
wind  tunnel  for  all  configurations.  There  is  a  increase  in  scatter  and  Cramer-Rao  bound 
magnitude  with  increasing  flap  deflection. 


5.5  Yaw  Rate  Derivatives 

The  yaw  rate  derivatives,  Cn^  and  Ci^,  are  plotted  against  a  in  figures  11  to  13.  The  AM- 
RL  low-speed  wind  tunnel  has  no  facility  to  measure  dynamic  derivatives,  and  empirical 
data  was  not  available  for  the  Pilatus  PC  9/A  in  approach  or  departure  configurations. 
As  mentioned  previously,  Cy^  was  constrained  to  its  a  priovi  value  for  the  maximum  like¬ 
lihood  analysis,  and  not  estimated  in  the  stepwise  regression  analysis.  Investigations  have 
previously  been  undertaken  [2,  Appendix  B]  to  assess  the  relative  contribution  of  each 
derivative  to  the  aircraft’s  forces  and  moments  to  ensure  that  those  set  to  their  a  priori 
values  were  small  contributors.  This  analysis  showed  the  yaw  rate  contribution  to  side 
force  to  be  small. 

Flight  test  estimates  of  yawing  moment  due  to  yaw  rate  derivative,  Cn^ ,  show  good  agree¬ 
ment  between  the  maximum  likelihood  and  stepwise  regression  techniques  for  all  config¬ 
urations.  Further,  the  small  Cramer-Rao  bounds  indicate  a  high  degree  of  confidence  in 
the  estimate,  despite  the  lack  of  corroborating  data.  The  figures  show  a  general  decrease 
in  the  average  magnitude  of  as  flap  deflection  is  increased. 

The  rolling  moment  due  to  yaw  rate,  Ci^,  whilst  being  small,  is  an  important  lateral- 
directional  cross-coupling  derivative.  Flight  test  estimates  of  Ci^  show  the  same  trends  as 
Cnr',  small  Cramer-Rao  bounds  indicating  a  high  level  of  confidence  in  the  estimate,  and 
a  general  decrease  in  the  magnitude  of  the  estimate  with  increasing  flap  deflection. 
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5.6  Roll  Rate  Derivatives 

Figures  14  to  16  show  the  roll  rate  derivatives,  Cnj,  and  Ci^,  plotted  against  a.  The 
side  force  due  to  roll  rate,  is  a  small  derivative  and  was  not  estimated  in  stepwise 
regression  or  maximum  likelihood  analysis. 

The  yawing  moment  due  to  roll  rate  derivative,  Cn^,  is  an  important  lateral-directional 
cross-coupling  term.  Maximum  likelihood  and  stepwise  regression  estimates  show  good 
agreement  with  low  Cramer-Rao  bounds  and  little  scatter.  The  average  magnitude  of 
decreases  as  flap  deflection  increases. 

The  roll  damping  derivative,  Ci^ ,  estimates  from  maximum  likelihood  and  stepwise  regres¬ 
sion  show  good  agreement,  with  small  Cramer-Rao  bounds.  The  magnitude  of  Ci^  does 
not  change  significantly  with  configuration. 


5.7  Lateral  Control  Derivatives 

Figures  17  to  19  show  the  lateral  control  derivatives,  and  plotted  against 

a.  The  yawing  moment  due  to  aileron  derivative,  ,  and  rolling  moment  due  to  rudder 
derivative,  ,  are  both  small  derivatives  which  were  not  estimated  in  either  stepwise 
regression  or  maximum  likelihood  analyses. 

Flight  test  estimates  of  show  a  higher  degree  of  scatter  than  is  present  in  estimates  of 
any  other  derivative,  but  are  centred  on  the  wind  tunnel  estimate.  The  maximum  likeli¬ 
hood  estimates  have  sizeable  Cramer-Rao  bounds,  indicating  a  reduced  level  of  confidence 
in  the  estimates,  while  the  stepwise  regression  estimates  have  quite  small  error  bands. 
There  is  a  general  increase  in  scatter  and  Cramer-Rao  bound  magnitude  with  increasing 
flap  deflection. 

Estimates  of  the  rudder  effectiveness  derivative,  >  show  small  Cramer-Rao  bounds  and 
error  bands,  indicating  a  high  level  of  confidence  in  the  estimates,  and  little  scatter.  Both 
majdmum  likelihood  and  stepwise  regression  estimates  are  larger  in  magnitude  than  the 
wind  tunnel  estimates.  There  is  no  significant  change  in  the  estimated  value  of  with 
aircraft  configuration. 

Estimates  of  the  aileron  effectiveness  derivative,  ,  show  little  scatter,  and  only  a  small 
increase  in  Cramer-Rao  bounds  with  increasing  flap  deflection.  Both  maximum  likelihood 
and  stepwise  regression  techniques  provide  estimates  that  are  consistently  larger  in  mag¬ 
nitude  than  the  wind  tunnel  estimates,  this  difference  tending  to  increase  with  increasing 
flap  deflection. 


5.8  Variation  of  Derivatives  with  Airbrake  Deployment 

The  effect  of  deployment  of  the  aircraft’s  airbrake  on  the  aerodynamic  characteristics 
of  the  aircraft  was  assessed  by  conducting  longitudinal  and  lateral  doublet  manoeuvres 
with  the  aircraft  in  clean  configuration  and  airbrake  extended.  Only  eight  flight  test 
manoeuvres  were  performed,  and  therefore,  the  results  from  maximum  likelihood  and 
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stepwise  regression  techniques  are  presented  as  averaged  values.  These  are  compared  with 
average  values  for  the  clean  aircraft  in  tables  3  to  4. 


Derivative 

Airbrake 

(Maximum  Likelihood) 

Airbrake 

(Stepwise) 

Clean 

Aircraft 

CNa 

9.23E-2 

8.59E-2 

8.53E-2 

Crria 

-6.28E-3 

-4.92E-3 

-6.53E-3 

Cm, 

-1.57E+1 

-1.95E+1 

-1.53E+1 

1.19E-3 

-4.67E-3 

1.42E-4 

Cms, 

-2.18E-2 

-2.02E-2 

-1.93E-2 

Table  3:  Comparison  of  averaged  longitudinal  flight  test  derivatives  with  airbrake  deployed, 
using  maximum  likelihood  and  stepwise  regression  techniques,  to  averaged  clean  aircraft 
results  [1]. 


Derivative 

Airbrake 

(Maximum  Likelihood) 

Airbrake 

(Stepwise) 

Clean 

Aircraft 

Cy, 

-1.65E-2 

-1.52E-2 

-1.32E-2 

Cn, 

1.68E-3 

1.94E-3 

1.44E-3 

Ch 

-1.45E-3 

-1.30E-3 

-1.40E-3 

Cnr 

-1.68E-1 

-2.03E-1 

-2.11E-1 

Ci^ 

1.32E-1 

1.24E-1 

1.37E-1 

Cnp 

-3.46E-2 

-2.36E-2 

-5.56E-2 

Cip 

-4.85E-1 

-4.45E-1 

-4.75E-1 

Cy^T 

2.03E-3 

1.72E-3 

3.38E-3 

Cns, 

-1.81E-3 

-1.90E-3 

-1.90E-3 

Cl,. 

-1.69E-3 

-1.68E-3 

-1.73E-3 

Table  4:  Comparison  of  averaged  lateral  flight  test  derivatives  with  airbrake  deployed, 
using  maximum  likelihood  and  stepwise  regression  techniques,  to  averaged  clean  aircraft 
results  [2]. 

The  deployment  of  the  airbrake  could  be  expected  to  have  its  most  significant  impact  in 
terms  of  an  incremental  addition  to  the  pitching  moment  and  axial  force  coefficients,  with 
little  impact  on  the  associated  derivatives.  However,  inspection  of  tables  3  and  4  show 
some  interesting  trends.  Table  3  shows  an  increase  in  magnitude  of  ,  Cmq  > 

and  a  decrease  in  magnitude  of  Cmai  with  airbrake  deployment.  The  other  longitudinal 
derivative,  is  not  consistently  affected  by  airbrake  deployment.  Table  4  shows  an 

increase  in  magnitude  of  both  Cy^  and  Cn^ ;  a  decrease  in  magnitude  of  the  yaw  damping 
term  Cur ;  and  a  significant  decrease  in  magnitude  of  Cup  and  Cy^^ . 

The  changes  in  the  lateral  stability  and  control  derivatives  are  particularly  interesting,  in 
that  airbrake  deployment  would  not  intuitively  be  expected  to  have  a  significant  impact 
on  them.  The  increase  in  magnitude  of  Cy^  is  a  result  of  the  increase  in  side  area  of 
the  aircraft  when  it  is  yawed  with  the  airbrake  deployed.  The  increase  in  magnitude  of 
the  directional  stability  term,  Cn^,  then  follows  because  the  airbrake  is  located  aft  of 
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the  centre-of-gravity,  and  therefore,  the  side  force  will  cause  a  restorative  moment.  The 
decrease  in  yaw  damping,  Cnr  >  can  not  be  attributed  directly  to  the  forces  and  moments 
generated  by  the  airbrake,  but  rather  to  its  influence  on  the  flow  over  the  aircraft’s  ventral 
fin.  Deployment  of  the  airbrake  will  significantly  disturb  the  airflow  downstream.  Hence 
the  ventral  fin  , which  is  located  behind  the  airbrake  on  the  underside  of  the  fuselage,  will 
experience  reduced  authority.  Whilst  there  is  a  significant  decrease  in  the  magnitude  of 
Cjip  this  derivative  has  only  a  small  impact  on  the  overall  handling  of  the  aircraft  and 
hence  any  changes  will  be  relatively  insignificant.  The  most  interesting  lateral  change 
due  to  airbrake  deployment  is  the  reduction  in  side  force  due  to  rudder  derivative,  , 
and  the  fact  that  this  is  not  accompanied  by  a  similar  decrease  in  the  rudder  control 
power,  ■  This  result  suggests  that  while  the  side  force  for  a  given  rudder  deflection  is 
reduced,  the  moment  it  produces  about  the  aircraft  centre-of-gravity  is  not,  a  result  that 
is  not  intuitive. 


5.9  Configuration  and  Altitude  Dependency 

Figures  2  to  19  show  graphically,  and  tables  6  to  9  summarise,  the  longitudinal  and  lateral 
derivatives  estimated  from  flight  test  for  the  Pilatus  PC  9/A  in  three  configurations:  gear 
deployed,  gear  deployed  with  take-off  flap,  and  gear  deployed  with  landing  flap. 

Figures  2  to  4  show  an  increase  in  the  Cramer-Rao  bound  magnitude  with  increasing 
flap  deflection,  particularly  in  plots  of  Cn^.  An  analysis  of  the  signal  variance  for  the 
a,  az,  P,  and  ay  flight  test  data  channels  was  performed  by  filtering  each  data  channel 
using  the  ‘filtfilt’  function  in  MATLAB©.  The  variance  of  the  original  data  channels  was 
then  calculated  by  comparison  to  the  filtered  data.  The  results  of  this  analysis,  averaged 
over  all  manoeuvres,  are  summarised  in  table  5  as  ratios  of  the  signal  ^variance  for  each 
configuration  to  the  signal  variance  for  the  gear  down  configuration  (^?— )■ 


a 

az 

ay 

Gear  down 

1.00 

1.00 

1.00 

1.00 

Take-off  flap 

1.37 

1.06 

2.15 

2.84 

Landing  flap 

3.49 

3.71 

7.54 

9.05 

Table  5:  Signal  variance  ratios  for  gear  down,  take-off  flap,  and  landing  flap  configurations. 

Inspection  of  the  table  shows  that  the  variance  ratio  of  the  flight  test  data  increases  with 
increasing  flap  deflection,  up  to  a  maximum  factor  of  9  for  ay  with  landing  flap.  This 
increase  in  signal  variance  with  flap  deflection  can  be  attributed  to  the  design  of  the  flaps 
on  the  PC  9/A.  The  aircraft  is  fitted  with  split  trailing  edge  flaps,  such  that  as  flap 
deflection  increases,  large  areas  of  separated  flow  are  produced  at  the  wing  trailing  edge. 
This  turbulent  separated  flow  buffets  the  empennage  of  the  aircraft  causing  the  whole 
airframe  to  vibrate.  It  is  these  vibrations  that  result  in  the  variance  shown  in  the  flight 
test  measured  data.  The  increase  in  variance  of  the  flight  test  data  manifests  itself  as  an 
increase  in  the  Cramer-Rao  bounds  of  the  flight  test  estimates  in  two  ways.  In  the  case  of 
the  a  and  fl  signals,  the  data  are  used  directly  by  the  maximum  likelihood  technique  and 
as  such  will  influence  the  quality  of  the  model  response.  For  the  acceleration  signals,  the 
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data  are  used  for  comparison  with  the  calculated  model  response,  and  as  such  any  variance 
will  manifest  itself  as  an  increase  in  error  between  the  model  and  the  flight  test  and  hence 
an  increased  Cramer-Rao  bound.  The  stepwise  regression  technique  is  not  affected  in  the 
same  way  by  signal  variance. 

Plots  of  Cy^  in  figures  8  to  10  show  a  divergence  between  flight  test  and  wind  tunnel 
estimates  as  flap  deflection  is  increased.  The  trend,  which  is  present  in  some  other  deriva¬ 
tives  to  a  lesser  extent,  could  be  the  result  of  several  factors.  As  described  previously,  the 
increased  flow  separation,  and  associated  impacts  on  the  flight  data,  that  occur  when  flaps 
are  deflected  could  lead  to  errors  in  the  flight  test  estimated  value  of  derivatives. 

Tables  6  to  9  show  the  mean  values  across  the  linear  angle-of-attack  range  (0°  -  6°)  of  the 
aerodynamic  derivatives  estimated  by  both  maximum  likelihood  and  stepwise  regression 
techniques,  for  each  configuration  and  altitude  band. 


Gear 
Down 
5000  ft 

Gear 

Down 

15  000  ft 

Take-off 

Flap 

5  000  ft 

Take-off 
Flap 
15000  ft 

Landing 

Flap 

5  000  ft 

Landing 

Flap 

15  000  ft 

8.83E-2 

8.62E-2 

9.51E-2 

9.16E-2 

9.40E-2 

9.35E-2 

Cma 

-1.03E-2 

-9.08E-3 

-1.15E-2 

-l.OlE-2 

-l.OOE-2 

-9.41E-3 

Cm, 

-1.54E+1 

-1.61E+1 

-1.14E-M 

-1.29E+1 

-1.49E+1 

-1.68E-1-1 

Cns^ 

2.62E-3 

1.62E-3 

2.99E-3 

2.60E-3 

1.37E-3 

1.63E-3 

n 

'-'ms. 

-2.34E-2 

-2.31E-2 

-2.32E-2 

-2.42E-2 

-2.36E-2 

-2.51E-2 

Table  6:  Summary  of  longitudinal  derivatives  estimated  using  maximum  likelihood. 


Gear 

Down 
5000  ft 

Gear 

Down 

15  000  ft 

Take-off 

Flap 

5  000  ft 

Take-off 
Flap 
15000  ft 

Landing 

Flap 

5  000  ft 

Landing 

Flap 

15  000  ft 

Cnc 

8.32E-2 

8.52E-2 

9.07E-2 

8.60E-2 

8.63E-2 

8.26E-2 

^rriot 

-8.75E-3 

-7.92E-3 

-8.80E-3 

-8.42E-3 

-7.87E-3 

-7.92E-3 

Cm, 

-1.74E+1 

-1.62E+1 

-1.58E+1 

-1.70E+1 

-1.77E+1 

-1.88E+1 

Cns^ 

-4.23E-3 

-4.23E-3 

-5.15E-3 

-2.75E-3 

-4.50E-3 

-3.55E-3 

Cms, 

-1.90E-2 

-1.78E-2 

-2.00E-2 

-2.12E-2 

-1.99E-2 

-2.06E-2 

Table  7:  Summary  of  longitudinal  derivatives  estimated  by  stepwise  regression. 

The  flight  test  estimates  of  the  longitudinal  derivatives  in  general  show  no  significant 
trend  with  configuration,  with  the  exception  that  the  maximum  likelihood  estimates  of 
and  Cms  show  a  decrease  and  increase  in  magnitude  respectively  with  increasing  flap 
deflection.  However,  many  of  the  lateral  derivative  flight  test  estimates  show  definite  trends 
with  flap  deflection.  Estimates  of  Cy^,  Cn^,  and  Cy^^  from  both  maximum  likelihood 
and  stepwise  regression  show  an  increase  in  magnitude  with  increasing  flap  deflection. 
Estimates  of  Cnj,,  and  (7/^  show  a  decrease  in  magnitude  with  flap  deflection.  The 
other  lateral  derivative  estimates  are  not  consistently  and  significantly  affected  by  flap 
deflection.  Figures  2  to  19  do  not  show  any  significant  altitude  dependency. 
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Gear 
Down 
5000  ft 

Gear 

Down 

15  000  ft 

Take-off 

Flap 

5  000  ft 

Take-off 
Flap 
15000  ft 

Landing 

Flap 

5  000  ft 

Landing 

Flap 

15  000  ft 

-1.75E-2 

-1.73E-2 

-2.05E-2 

-1.96-2 

-2.29E-2 

-2.22E-2 

Cn, 

1.59E-3 

1.54E-3 

2.01E-3 

1.89E-3 

2.18E-3 

2.22E-3 

-1.60E-3 

-1.50E-3 

-1.43E-3 

-1.62E-3 

-1.17E-3 

-1.48E-3 

Cnr 

-1.97E-1 

-1.99E-1 

-1.52E-1 

-1.48E-1 

-1.39E-1 

-9.57E-2 

Cl. 

1.91E-1 

1.86E-1 

5.79E-2 

5.91E-2 

1.05E-1 

5.54E-2 

Cn, 

-7.41E-2 

-7.11E-2 

-4.34E-2 

-3.58E-2 

-4.27E-2 

-2.88E-2 

Cl, 

-5.27E-1 

-5.00E-1 

-4.89E-1 

-5.63E-1 

-4.54E-1 

-5.63E-1 

Cys^ 

4.17E-3 

4.41E-3 

4.72E-3 

4.97E-3 

5.23E-3 

7.54E-3 

Cns^ 

-2.12E-3 

-2.05E-3 

-2.01E-3 

-1.82E-3 

-1.99E-3 

-1.98E-3 

-1.90E-3 

-1.78E-3 

-1.94E-3 

-2.04E-3 

-1.90E-3 

-2.05E-3 

Table  8:  Summary  of  lateral  derivatives  estimated  using  maximum  likelihood. 


Gear 
Down 
5000  ft 

Gear 

Down 

15  000  ft 

Take-off 

Flap 

5  000  ft 

Take-off 
Flap 
15000  ft 

Landing 

Flap 

5  000  ft 

Landing 

Flap 

15  000  ft 

Cyd 

-1.59E-2 

-1.67E-2 

-1.81E-2 

-1.80E-2 

-2.14E-2 

-2.09E-2 

Cn, 

1.67E-3 

1.61E-3 

2.09E-3 

2.05E-3 

2.32E-3 

2.35E-3 

Cl, 

-1.58E-3 

-1.48E-3 

-1.34E-3 

-1.37E-3 

-1.15E-3 

-1.30E-3 

Cnr 

-1.96E-1 

-2.07E-1 

-1.60E-1 

-1.67E-1 

-1.54E-1 

-1.22E-1 

Clr 

1.20E-1 

1.34E-1 

1.17E-1 

1.48E-1 

1.27E-1 

1.36E-1 

Cn, 

-5.98E-2 

-5.85E-2 

-2.47E-2 

-3.72E-2 

-2.37E-2 

-2.22E-2 

Cl, 

-5.02E-1 

-4.69E-1 

-4.65E-1 

-4.56E-1 

-4.56E-1 

-4.59E-1 

CYs, 

2.76E-3 

2.78E-3 

2.79E-3 

2.99E-3 

2.96E-3 

2.92E-3 

Cns^ 

-2.14E-3 

-2.09E-3 

-2.04E-3 

-2.10E-3 

-2.07E-3 

-2.06E-3 

Ci^ 

-1.90E-3 

-1.78E-3 

-1.92E-3 

-1.89E-3 

-1.92E-3 

-1.95E-3 

Table  9:  Summary  of  lateral  derivatives  estimated  by  stepwise  regression. 

6  Conclusions 

The  longitudinal  and  lateral  stability  and  control  derivatives  of  the  PC  9/ A  in  approach 
and  departure  configurations  have  been  determined  from  flight  test  data  measurements 
using  maximum  likelihood  and  stepwise  regression  estimation  techniques.  Where  available, 
comparisons  have  been  made  with  derivative  estimates  from  wind  tunnel  testing. 

Flight  test  estimates  of  the  normal  force  due  to  angle-of-attack  derivative,  Cno,  showed 
an  increase  in  magnitude  with  increasing  flap  deflection,  and  good  agreement  with  wind 
tunnel  results  for  all  configurations.  Pitch  stifiness  derivative,  Cmai  estimates  from  flight 
test  were  consistently  smaller  in  magnitude  than  wind  tunnel  results.  However,  the  same 
trend,  of  increasing  pitch  stiffness  with  flap  deflection,  was  demonstrated  in  flight  test  and 
wind  tunnel.  Flight  test  estimates  of  are  scattered  about  zero,  but  in  all  cases  are 
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less  than  the  wind  tunnel  result.  Maximum  likelihood  estimates  of  elevator  effectiveness, 
Cms  )  match  well  with  wind  tunnel  results  for  all  configurations,  while  stepwise  regression 
estimates  are  consistently  smaller  in  magnitude.  No  wind  tunnel  or  empirical  results 
for  pitch  stiffness  were  available  for  comparison  with  flight  test  estimates.  However,  the 
flight  test  estimates  show  a  general  trend  of  decreasing  pitch  stiffness  with  increasing  flap 
deflection. 

Flight  test  estimates  of  angle-of-sideslip  derivatives,  Cy^  and  ,  and  rudder  derivative, 
Cy^  ,  show  good  agreement  with  wind  tunnel  estimates  for  all  conflgurations.  Flight  test 
estimates  of  directional  stability  derivative,  Cn^,  are  consistently  smaller  in  magnitude 
than  the  wind  tunnel  results,  for  all  conflgurations.  Flight  test  estimates  of  derivatives  0^^ 
and  Ci^  are  greater  in  magnitude  than  wind  tunnel  results.  No  wind  tunnel  or  empirical 
data  were  available  for  comparison  of  lateral  dynamic  derivatives.  However,  rolling  moment 
derivatives,  Ci^  and  Ci^  were  largely  unaffected  by  configuration  changes,  while  yawing 
moment  derivatives,  Cn^  and  Cn,,  demonstrated  a  trend  of  decreasing  magnitude  with 
increasing  flap  deflection. 

In  general,  the  flight  test  results  showed  no  significant  trend  with  altitude.  However,  all 
flight  test  estimated  derivatives  showed  a  trend  of  increased  scatter  and  larger  Cramer-Rao 
bounds  with  increasing  flap  deflection.  Analysis  of  flight  test  data  showed  that  this  trend 
can  be  attributed  to  changes  in  signal  variance  due  to  increased  airframe  bufi'et  with  flap 
deployment. 

Comparisons  of  flight  test  estimates  of  stability  and  control  derivatives,  with  airbrake 
deployed,  to  those  of  the  clean  aircraft  showed  significant  changes  both  longitudinally  and 
laterally.  However,  the  limited  number  of  test  manoeuvres  performed  in  this  configuration 
places  a  caveat  on  any  conclusions  drawn  about  the  influence  of  the  airbrake. 
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Figure  1:  Air  data  boom  installation. 
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Figure  2:  PC  9/A  angle-of-attack  and  pitch  rate  derivatives,  gear  down. 
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Figure  4-'  PC  9/A  angle- of- attack  and  pitch  rate  derivatives,  landing  flap. 
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Figure  5:  PC  9/A  longitudinal  control  derivatives,  gear  down. 
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Figure  6:  PC  9/ A  longitudinal  control  derivatives,  take-off  flap. 
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Figure  7:  PC  9/A  longitudinal  control  derivatives,  landing  flap. 
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Figure  8:  PC  9/A  angle-of -sideslip  derivatives;  gear  down. 
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Figure  10:  PC  9/ A  angle- of- sideslip  derivatives,  landing  flap. 
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Figure  12:  PC  9/A  yaw  rate  derivatives,  take- ojf  flap. 
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Figure  14-  PC  9/ A  roll  rate  derivatives,  gear  down. 
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Figure  15:  PC  9/ A  roll  rate  derivatives,  take-off  flap. 


(per  radian)  C  ^  (P^r  radian) 


DSTO-TR-0479 


O  Flight  Test  (pEst  estimates)  5000  ft 
□  Flight  Test  (swr  estimates)  5000  ft 
•  Flight  Test  (pEst  estimates)  15  000  ft 
■  Flight  Test  (swr  estimates)  15  000  ft 


Figure  16:  PC  9/ A  roll  rate  derivatives,  landing  flap. 
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Figure  1 7:  PC  9/A  lateral  control  derivatives,  gear  down. 
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Figure  18:  PC  9/A  lateral  control  derivatives,  take-off  flap. 
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Figure  19:  PC  9/ A  lateral  control  derivatives,  landing  flap. 
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Appendix  A:  Mass,  Cent  re- of- Gravity  and 
Mass  Moments-of-Inertia 


Equations  for  the  mass,  centre-of-gravity,  and  mass  moments-of- inertia  of  the  flight  test  ar¬ 
ticle  in  clean  configuration  were  derived  as  functions  of  fuel  mass  in  reference  [1,  Appendix 
A].  These  equations  are  summarised  below. 


M  =  1866.1  +  Mf  (kg) 

(Al) 

^  8015.5-I-4.178MP 

•^c.g.  —  ism.l+MF 

(A2) 

^XX 

=  2505.9  + 6.177  (M 

-  1866.1) 

(kg.m^) 

(A3) 

Iyy 

=  6622.2  +  0.033  (M 

-  1866.1) 

(kg.m^) 

(A4) 

Izz 

=  8467.1  +  6.188  (M 

-  1866.1) 

(kg.m^) 

(A5) 

IxY 

=  49.0 +  0.0057  (M- 

-  1866.1) 

(kg.m^) 

(A6) 

Ixz 

=  196.9  +  0.0041  {M 

-  1866.1) 

(kg.m^) 

(A7) 

Iyz 

=  3.0  + 0.0007  (M  - 

■  1866.1) 

(kg.m^) 

(A8) 

During  the  flight  test  program,  the  clean  aircraft  longitudinal  centre-of-gravity  position 
varied  between  4.27m  (24.2%MAC)  and  4.30m  (25.6%MAC).  The  lateral  and  vertical 
centre-of-gravity  positions  were  assumed  to  be  invariant  with  fuel  usage,  and  were  fixed 
at  values  of  0.024m  and  -2.2m,  respectively,  relative  to  the  aircraft  datum. 

For  the  purposes  of  analysing  the  approach  and  departure  configuration  flight  test  data, 
the  above  equations  for  centre-of-gravity  location  and  mass  moments-of-inertia  have  been 
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modified  to  account  for  the  effects  of  deployment  of  the  landing  gear.  The  effect  on 
centre-of-gravity  location  resulting  from  fiap  and  speed  brake  deployment  was  assumed  to 
be  negligible. 

Tables  Al,  A2,  and  A3  detail  the  clean  aircraft  basic  mass  distribution,  retracted  landing 
gear  mass  distribution,  and  deployed  landing  gear  mass  distribution,  respectively.  This 
data  is  sourced  from  references  [3,  4,  18,  19].  All  distances  are  relative  to  the  aircraft 
datumn  located  3  metres  forward  of  the  engine  firewall  and  2  metres  below  the  forward 
fuselage  reference  line.  The  landing  gear  has  been  modelled  as  point  masses  located  at  the 
end  of  each  leg. 


Mass  (kg) 

Moment  Arm  (m) 

Basic  Aircraft 

1784.5 

4.306 

Pilot 

81.6 

4.061 

Fuel 

Variable 

4.178 

Table  Al:  Flight  test  aircraft  mass  distribution  [3,  41- 


Mass  (kg) 

X  (m) 

y  (m) 

^  (m) 

Nose  gear 

42.0 

3.45 

0.00 

1.90 

Starboard  main  gear 

45.3 

4.62 

0.50 

1.90 

Port  main  gear 

45.3 

4.62 

-0.50 

1.90 

Table  A2:  Retracted  landing  gear  mass  and  moment  arm. 


Mass  (kg) 

X  (m) 

y  (m) 

^  (m) 

Nose  gear 

42.0 

2.28 

0.00 

1.08 

Starboard  main  gear 

45.3 

4.62 

1.27 

1.23 

Port  main  gear 

45.3 

4.62 

-1.27 

1.23 

Table  A3:  Deployed  landing  gear  mass  and  moment  arm. 

The  PC  9/A  main  gear  deploys  and  retracts  laterally  with  little  or  no  longitudinal  motion. 
Therefore,  the  effect  of  its  deployment  on  the  longitudinal  centre-of-gravity  location  of  the 
aircraft  was  assumed  to  be  negligible.  However,  the  effect  of  nose  gear  deployment  on 
longitudinal  centre-of-gravity  location  was  taken  into  account.  The  longitudinal  location 
of  the  aircraft  centre-of-gravity,  with  landing  gear  deployed,  is  given  by  equation  A9. 


m. 


X, 


ng 


C-g-depto 


^c.g. 


f^tot 


^C.g.ng^^t  ^C.g.ng^^p 


)■ 


(A9) 


Where  is  the  mass  of  the  nose  gear,  Xc.g.^g^^^  is  the  centre-of-gravity  of  the  nose 
gear  in  the  retracted  position,  is  the  centre-of-gravity  of  the  nose  gear  in  the 

deployed  position,  and  x^.g.  is  the  clean  aircraft  centre-of-gravity  location.  Substitution  of 
the  appropriate  data  into  equation  A9  yields, 
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^C.g.deptot  ■“  ^<=.3. 


48.972 

mtot 


(m) 


(AlO) 


This  represents  a  forward  longitudinal  c.g.  shift  of  approximately  0.5%c. 

Similar  equations  were  developed  to  determine  the  change  in  the  vertical  c.g.  position. 


^C-.g-deptot 


c.g. 


mtot  [  c.g. 


™3dep 


•'C.g.mgret 


)  ■  (All) 


As  the  vertical  centre-of-gravity  location  of  the  clean  aircraft  was  assumed  constant,  all 
terms  except  mtoti  iii  equation  All,  were  known.  Substitution  of  the  relevant  data  yields, 


=  2.200 


107.868 

mtot 


(A12) 


This  represents  a  vertical  c.g.  shift  of  approximately  2.2%c  downwards. 

To  evaluate  the  effect  of  landing  gear  deployment  on  the  mass  moments-of-inertia  of  the 
aircraft,  the  data  from  tables  A2  and  A3  were  used  to  evaluated  the  inertia  of  the  individual 
gear  components  in  both  the  retracted  and  deployed  positions.  The  difference  between 
the  deployed  values  and  the  retracted  values  represents  the  increment  that  was  added  to 
the  clean  aircraft  values  yielding  inertial  data  for  the  aircraft  in  approach  and  departure 
configurations.  These  increments  are  detailed  in  table  A4. 


Parameter 

kg.m^ 

AJxx 

249.70 

AIyy 

266.70 

A/zz 

263.80 

AIxy 

1.16 

^Ixz 

64.95 

AIyz 

2.70 

Table  AJti  Mass  moment- oj-inertia  increments  due  to  landing  gear  deployment. 

The  inertia  increments  due  to  gear  deployment  were  added  to  equations  A3  to  A8.  This 
yielded  the  gear  deployed  inertia  equations. 


Ixx, 


dePtot 


^yYdep,,t 


=  2505.9  +  6.177  (M-  1866.1) 


=  6622.2  +  0.033  (M  -  1866.1) 


(kg.m^) 

(A13) 

(kg.m^) 

(A14) 
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=  8467.1  + 6.188  (M-  1866.1)  (kg.m^)  (A15) 

=  49.0  + 0.0057  (M-  1866.1)  (kg.m^)  (A16) 

=  196.9  +  0.0041  (M  -  1866.1)  (kg.m^)  (A17) 

=  3.0  + 0.0007  (M-  1866.1)  (kg.m^)  (A18) 
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